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ABSTRACT 


QCD predicts a phase transition between hadronic matter and a Quark 
Gluon Plasma at high energy density. The Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory is a new facility dedicated to 
the experimental study of matter under extreme conditions. I will discuss 
the first round of experimental results from colliding heavy nuclei at RHIC 
and our current understanding of the state of matter generated in such col¬ 
lisions, concentrating on partonic energy loss and jet quenching as a probe 
of the medium. 


*pmjacobs@lbl.gov 


1 



1 Introduction 


At high temperature or baryon density, hadronic matter dissolves into a soup of its con¬ 
stituent quarks and gluons. For an asymptotically free field theory such as QCD, the 
state of matter at high energy density is simple!: long range (low momentum) inter¬ 
actions are screened, and short range (high momentum) interactions are weak, lead¬ 
ing to an ideal gas equation of state in the high energy density limit. At temperature 
T Aqcd matter is a gas of deconfined, weakly interacting quarks and gluons (the 
“Quark-Gluon Plasma”, or QGP), whereas at T Aq C d quarks and gluons are con¬ 
fined and matter consists of strongly interacting hadrons. 

The QCD phase diagram has a 
complex structure (Fig. [1]). At 
low temperature and baryon density the 
phase is hadronic (confined phase) and 
chiral symmetry is broken. Color¬ 
superconducting and other phases may 
exist at high baryon density and low tem¬ 
perature,! whereas at high temperature 
the quarks and gluons are deconfined and 
chiral symmetry is restored. The early 
universe descended from high T at ex¬ 
tremely small fiB- Neutron star cores 
have high /i B and very low T. 

First-principles calculations of finite 
temperature QCD can only be carried out numerically on the lattice.! Fig. |2j shows 
a recent lattice calculation of the energy density e as function of temperature for two- 
and three-flavor QCD. e exhibits a sharp rise in the vicinity of the critical temperature 
T c , indicating a rapid change in the density of underlying degrees of freedom. The ideal 
gas Stefan-Boltzmann limit e S B has not yet been achieved at T ~ 4Tc. Putting in phys¬ 
ical values, Tc ~ 175 MeV, resulting in critical energy density ec = (6 ± 2)T C 4 ~ 1 
GeV/frn 3 . This value should be kept in mind for comparison to conditions achieved in 
laboratory experiments. 

The order of the deconfinement phase transition can be determined in some limiting 
cases.! It is first order for pure gauge and for three light quarks, second order for two 
light and one heavy quark. For physical quark masses the order of the transition, or 



Fig. 1. Example of the QCD phase diagram 
(temperature T vs. baryochemical poten¬ 
tial // rj ) for two massless flavors.! Chiral 
symmetry is broken in the hadronic phase. 
“2SC” is a color-superconducting phase. 
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indeed whether it is a smooth cross over, has not been determined. The extension of 
lattice calculations to /i B 7 ^ 0 is a long-standing problem, but there has been significant 
recent progress in determining the phase boundary and equation of state for finite // B .0 

In the early universe, the confinement 
transition (QGP—>hadrons) at very low 
/iB occurred about 10 microseconds after 
the big bang.i A strongly first order phase 
transition may have generated primordial 
black holes, strange quark nuggets, or lo¬ 
cal baryon asymmetries affecting primor¬ 
dial nucleosynthesis, though no relics of 
this transition have been observed thus 

tt t t f j far. In the cores of neutron stars, a QGP 

Fig. 2. Lattice calculation of energy den- ^ 

.. , r , ,,, n phase at low temperature and high u B 

sity vs temperature for two- and three-flavor 1 1 fc ' 

| may generate observable millisecond pul¬ 

sar phenomena.! 

It is natural to ask whether the deconfinement and chiral symmetry restoration tran¬ 
sitions can be studied in accelerator-based experiments. We require a system having 
temperature of order the pion mass in equilibrium over a large volume. The best possi¬ 
bility to accomplish this is through the collision of heavy nuclei at the highest possible 
energy. 

However, nuclear collisions at high energy are highly dynamic. Even if a hot, equi¬ 
librated system is created early in the collision, it will immediately expand and cool. If 
a QGP is formed it will hadronize after a brief period, and signatures from the decon- 
fined phase may be masked by those from the hot, interacting hadron gas. The central 
issue is to find and measure those experimental observables that are sensitive to the 
state of matter early in the collision.! 

Broadly speaking, several epochs in the evolution of high energy nuclear collisions 
can be sketched. Immediately following the collision is a brief formation time r < 1 
fm/c, during which large momentum transfer (short-distance) processes occur. It is also 
during this period that the highest energy density is achieved and the QGP may result. 
After an expansion time of perhaps a few fm/c the plasma hadronizes into a dense, in¬ 
teracting hadron gas. Further expansion and cooling causes the temperature to fall to 
the point at which inelastic collisions among hadrons are no longer common {chemical 
freezeout) and the relative populations of the various long-lived hadron species are es- 
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tablished. Elastic collisions occur until kinetic freezeout, after which point the hadrons 
fly effectively undisturbed to the detectors. Thus, the relative populations of stable 
hadrons reflect the conditions at chemical freezeout, whereas their momentum spectra 
reflect the conditions at kinetic freezeout. More penetrating probes (dileptons, direct 
photons, jets, heavy quarks,...) may carry information on the conditions around the 
time of their formation earlier in the collision. 

Collisions of heavy nuclei at ultra-relativistic energies have been studied in fixed 
target experiments at the Brookhaven AGS and CERN SPS for the past 15 years. In 
1999, the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
was commissioned, bringing into operation the first facility largely dedicated to the 
study of matter at high energy density. 

I will first describe the RHIC machine and experiments and then discuss some of the 
main results from its heavy ion physics program. I will briefly touch on “soft physics” 
observables {px <~ 2 GeV/c) and summarize what has been learned from them. The 
bulk of the review concentrates on the major new development at RHIC: jet production 
and indications of partonic interactions with dense matter, which may directly probe 
the energy density achieved early in the collision. 

2 The Relativistic Heavy Ion Collider 

RHIC consists of two independent superconducting rings, 3.8 km in length. It has enor¬ 
mous flexibility in beam masses and energies, with capability to collide gold ions from 
a /s nn =20 to 200 GeV per nucleon pair, protons up to i/s=500 GeV, and asymmetric 
systems, most importantly protons and deuterons with heavy nuclei. The top CM en¬ 
ergy for heavy nuclei is a factor 10 larger than for the fixed target experiments, extend¬ 
ing significantly the statistical and transverse momentum reach of many observables 
and opening up new channels. RHIC is also the first polarized proton collider, opening 
new opportunities to study the spin content of the proton, in particular the contribution 
of the gluon at low ^Bj. More details of the machine and its first year performance can 
be found in Ref. (j9j). 

Design luminosity for Au+Au at 200 GeV is C = 2 • 10 26 cm _2 sec _1 , giving an 
interaction rate of about 2 kHz. While this luminosity appears to be tiny relative to 
other modem colliders, recall that the rate for hard processes such as jet production in 
nuclear collisions scales as ~ A 2 (A=atomic mass), so that hard process rates in Au+Au 
at design luminosity are the same as at a proton collider with C ~ 10 31 cm _2 sec _1 . 
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The design luminosity for p+p collisions at 500 GeV is 2 • 10 32 cm” 2 sec _1 , with an 
interaction rate of about 8 MHz. Design polarization for p+p is 70%. 

RHIC has six intersection regions, of which four are currently instrumented with ex¬ 
periments.0 PHENIX consists of an axial field magnet and four independent spectrom¬ 
eters: two at midrapidity containing tracking, ring imaging Cerenkov counters, time of 
flight, and electromagnetic calorimetry, which are optimized for precision lepton, pho¬ 
ton and hadron measurements, and two forward muon arms. STAR has conventional 
collider detector geometry, with a large solenoidal magnet, Time Projection Cham¬ 
ber for tracking, large coverage EM calorimetry, and an inner silicon-based tracker. 
STAR is designed for hadron, jet, lepton and photon measurements over large accep¬ 
tance, as well as studies of event-wise fluctuations in high multiplicity nuclear colli¬ 
sions. BRAHMS consists of two small acceptance spectrometers for inclusive identi¬ 
fied hadron measurements over wide phase space. PHOBOS has very wide phase space 
coverage for charged particles using silicon detectors, and a mid-rapidity spectrometer 
based on a dipole magnet, silicon tracking and Time of Flight. STAR and PHENIX 
each have about 450 collaborators, whereas PHOBOS and BRAHMS each have fewer 
than 100 . 

RHIC had a brief commissioning run in 1999. The data reported here are from 
a Au+Au run at v /snn= 130 GeV in 2000 and Au+Au and polarized p+p (~15% po¬ 
larization) runs at ^/s N n=200 GeV in 2001-2. The integrated luminosity is about 80 
/ib” 1 for Au+Au and 1 pb ” 1 for p+p. A comprehensive view of the physics program 
of RHIC and what has been achieved thus far can be found in the proceedings of the 
recent Quark Matter conferences. m 



vs BBC analog response 


Fig. 3. Event characterizatiorS correla¬ 
tion of ZDC signal (vertical) with multiplic¬ 
ity. 


Nuclei are extended objects, and high 
energy nuclear collisions can be charac¬ 
terized experimentally as head-on (“cen¬ 
tral”) or glancing (“peripheral”). For 
this purpose all experiments have iden¬ 
tical forward hadron calorimeters (Zero 
Degree Calorimeters, or ZDCs0), situ¬ 
ated ±18 m downstream of the inter¬ 
action region beyond the first machine 
dipole magnets. The ZDCs measure 
the forward-going “spectator” neutrons 
which did not scatter significantly. The 
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correlation of the ZDC signal with the event multiplicity serves to characterize the 
geometry of the event. Peripheral collisions have few spectator neutrons and low multi¬ 
plicity. More central collisions generate more spectator neutrons and higher multiplic¬ 
ity, while the most central, head-on collisions again generate few spectator neutrons 
and the highest multiplicity. Fig. |3] shows this correlation from the PHENIX exper¬ 
iment, the other experiments use similar distributions. Based on this correlation the 
total geometric cross section can be divided into bins containing a certain percentile of 
“centrality”, as shown by the lines on the figure. 

Related quantities which also characterize the centrality of the event are the num¬ 
ber of incoming participating nucleons (N part ), and the number of equivalent binary 
nucleon-nucleon collisions (N binary ). N part and N binary are calculated for each central¬ 
ity bin using the Glauber formalism and a realistic model of the nuclear geometry .111 
Phenomenologically it has been found that total particle production scales roughly as 
Npart, whereas the rate of hard processes will scale as N binary in the absence of nuclear 
effects. Such scaling rules can be used to uncover the new physics present in nuclear 
collisions. For instance, the violation of N binary scaling at high p T indicates significant 
effects of the nuclear medium on high p T processes. I will return to this point below. 


3 Soft Physics 

Fig. [Tj from the PHOBOS CollaborationS shows charged particle multiplicity distribu¬ 
tions dN ch /dr] vs //0 over the full RHIC phase space, for Au+Au collisions at all three 
collision energies studied so far. The longitudinal phase space growth with increasing 
energy is apparent, as well as the increase in multiplicity for more central collisions. 
For the most central collisions at x /s NN =130 GeV, 4200 charged particles are produced 
in the full phase space. 

Figure [4] exhibits a central plateau near y ~ 0, indicating that the system has ap¬ 
proximate longitudinal boost invariance. More detailed considerations of identified par¬ 
ticle spectra show that boost invariance holds over a rather smaller region Ay ~ 1.0 
Bjorken studied boost invariant hydrodynamics0 and derived a useful pocket formula 
for the energy density achieved in the central region: 


1 1 . .3 dN ch 

Bj itRa 2 t dy ttR a 2 t ^ T 2 dy 


( 1 ) 
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where Et and dN c h/drj are the transverse energy 
and multiplicity density (Fig. |]). R -\ is the nuclear 
radius and r is the formation time, typically taken as 
~ 1 fm/c, after which the hydrodynamic description 
is valid. PHENIX measured Et for central Au+Au 
collisions at v /skn= 130 GeV and derived an energy 
density of 4.6 GeV/fin 3 ,@ well above 1 GeV/fm 3 
required for the deconfinement transition as calcu¬ 
lated on the lattice. Conditions appear to be favor¬ 
able to form a QGP in nuclear collisions at RHIC, 
though under the assumption that an equilibrated 
system is established after a time ~ 1 fm/c. 

Scaling of global observables such as multiplic¬ 
ity with collision centrality and y/s has been dis¬ 
cussed extensively in the context of gluon saturation 
(e.g. Ref. |15] and references therein). A projectile 
parton at sufficiently low x interacts coherently with 
all partons in the target that lie within a transverse 
area 1 /Q 2 , with cross section a ~ a s /Q 2 . If the 
target is a heavy nucleus, its finite thickness will in¬ 
crease the parton density seen by the probe: nuclei are “blacker” than protons. The 
Q 2 scale at which the nucleus becomes dense to the probe is given by the saturation 
condition 1 


Fig. 4. dN ch /dr] for various col¬ 
lision centralities over the full 
RHIC phase space.li! 


Qt ~ 


2 xG a (x, Q 2 ^ 


7:R 2 a 


A 3 . 


( 2 ) 


where xGa ~ A is the nuclear gluon distribution and Ra ~ Az is the nuclear radius. 
For central Au+Au collisions at RHIC Q 2 ~ 2 GeV 2 , so that a s {Q 2 ) is small: in 
the saturation regime the coupling is weak while the gluon density is high, leading to 
classical, nonlinear dynamics. The growth of Q 2 with target thickness leads e.g. to 
definite predictions for the dependence of multiplicity on collision centrality .0 

Fig. |5(a)| from PHOBOS illustrates this dependence for Au+Au collisions at 130 
and 200 GeV (central collisions correspond to large number of participants N part , pe¬ 
ripheral to small N part ). The data are compared to both the saturation model and a more 
conventional “two component” model, which parameterizes the centrality dependence 
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(a) Multiplicity^ 
centrality Jiil 


0) per participant pair vs. 


(b) Universal scaling of total multiplicity vs. 
■s/s for e + e _ , p+p and Au+Au.@ 


Fig. 5. 


of multiplicity in terms of soft production (scaling with N part ) and minijet production 
(scaling with Nbmary)- The data are consistent with both models and cannot resolve the 
small difference between them. While this is a very active area of research, no clear 
evidence of saturation effects in nuclear collisions has been found yet at RHIC. 

Fig. |5(b)| from PHOBOS shows a universal scaling of the total multiplicity with 
a/s for e + e“, p+p and Au+Au.@ (The p+p data are corrected for energy carried by 
the leading particles.) This is at first sight surprising, given the very different nature of 
the various systems, but may point to universal features of particle production in high 
energy collisions (though it should also be noted that the evolution with yfs of the mean 
transverse momentum shows no such scaling among the various systemsEl). 

As argued above, the relative population of the stable hadrons is fixed at chemi¬ 
cal freezeout. Equilibrium at that point should be evident in the measured population 
ratios. Fig. |] shows population ratios measured at mid-rapidity for a wide variety of 
non-strange and strange hadrons, with values varying over three orders of magnitude. 
Also shown are results from a statistical model fit to the data.il The model consists of a 
partition function containing contributions from all baryons and mesons with mass less 
than 2 GeV, with two free parameters (chemical freezeout temperature T che m and bary- 
ochemical potential /i B ) and constraints due to strangeness and charge conservation. 
The resulting fit from this and similar models is excellent, with extracted parameters 
Tchem ~ 170 MeV and /x B ~ 25 — 50 MeV, close to the critical temperature for the 


8 
















1 


10 


10 



Fig. 6. Population ratios of stable hadrons and statistical model fit 


.Eill 


deconfinement phase transition calculated on the lattice. A hadronic gas cannot ex¬ 
ist at a higher temperature, suggesting that equilibrium is in fact established prior to 
hadronization, in the deconfined phase. 

I cannot hope to address all the soft physics observables that have been studied 
at RHIC. To conclude this section I will summarize a few of the main results, more 
detailed discussions can be found in Refs. (|12j,[r8[). 

• Low baryon density: antibaryon/baryon ratios at midrapidity are ~0.6-1.0. The 
system is close to baryon free, similar to the early universe, but not precisely so. 
Finite baryon number is transported Ay ~ 5.5 rapidity units from the beam. 

• There are strong indications that hadronic chemical equilibrium has been achieved 
at a temperature near the lattice critical temperature. 

• Hydrodynamic calculations describe well the transverse momentum spectra and 
“elliptic flow” (see below). The mass dependence of these observables, which is 
a sensitive test of the hydrodynamic picture, is described in detail. 

• The energy density achieved early in the collision is estimated to be ~ 5 GeV/fm 3 . 

• Identical two-particle correlations (Hanbury Brown-Twiss correlations) are sen¬ 
sitive to the space-time evolution of the source. The extracted source radii and 
duration of freezeout show no significant increase relative to lower energy nuclear 
collisions. Such correlations measure only a piece of a dynamically expanding 
source and these results may indicate a very explosive expansion resulting from 
high early pressure. 


9 




4 High p i Hadrons and Jets 


The increase in p snn for nuclear collisions at RHIC relative to fixed target experiments 
opens up new channels to probe the dense medium generated in the collision. Jets 
with Et ~ 40 GeV and higher are produced in sufficient numbers to provide robust 
observables. The measurement of jets in nuclear collisions poses a special problem, 
however: while the presence of a hard scattering in a nuclear collision can be detected 
(though in a biased way) via high px leading hadrons, the huge soft multiplicities seen 
in Fig. [4] contaminate any finite jet cone, spoiling the jet energy measurement. At 
sufficiently high E T this effect may be minor (e.g. nuclear collisions at the LHC), 
but for the jet E T currently accessible at RHIC it is fatal. We therefore restrict our 
considerations to leading particles and their correlations. We show below that hadrons 
with px >~ 4 GeV/c are produced dominantly from jet fragmentation, even in the most 
central Au+Au collisions. 

Twenty years ago BjorkcnEI proposed that hard scattered partons in nuclear col¬ 
lisions could provide a sensitive probe of the surrounding medium. The energy loss 
dE/dx due to elastic scattering of the partons in a Quark Gluon Plasma depends on the 
temperature as Tp lasma and results in a suppression of the observed rate of jets or their 
leading hadrons at fixed px- Later work showed that dE/dx from elastic scattering is 
negligible but that radiative energy loss in dense matter could be considcrablc.^Sill 
The energy loss is directly sensitive to the gluon density of the medium, p g i ue . While 
not a direct signature of deconfinement, measurement of p g i ue that is substantially larger 
than in cold nuclear matteril is incompatible with the presence of a hadronic medium, 
thus large energy loss serves as an indirect signature of deconfinement. 

There are currently three sets of measurements which address the question of par- 
tonic energy loss in dense matter, which I will discuss in turn: suppression of inclusive 
spectra, elliptic flow (azimuthal asymmetry) in non-central events, and correlations of 
high p T hadron pairs. 


4.1 High p T : Suppression of Inclusive Spectra 

In p+p collisions at RHIC energies, hadrons with p T ~ 4 GeV/c typically carry 75% 
of the energy of their parent jet, leading to the possibility that partonic energy loss in 
Au+Au collisions is reflected in the suppression of leading hadrons. Ell Hadron suppres¬ 
sion is measured via the nuclear modification factor: 
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Fig. 7. Inclusive charged hadron invariant distributions for various event centralitiesH 


Raa(pt ) 


da AA /dydp T 2 

( N binary ) d(Tpp / dy djfj'~ 


(3) 


where a pp /dyd 2 p T is the inclusive cross section measured in elementary nucleon-nucleon 
collisions and (N binary ) accounts for the geometric scaling from the elementary to nu¬ 
clear collision. Raa(pt ) is normalized to be unity if Au+Au collisions are an incoher¬ 
ent superposition of p+p collisions. In addition to partonic energy loss, Raa{pt ) may 
be altered by nuclear effects such as gluon shadowing, which will reduce R A a{pt), 
and the Cronin effect (multiple soft scattering in the initial state), which will increase 
it. These effects must be disentangled using measurements in simpler systems (in par¬ 
ticular p+Au). 

First indication of a strong suppression of high pt hadrons in nuclear collisions was 
reported from the 130 GeV data by both PHENIX and STAR.il The 200 GeV data have 
much higher statistics and both collaborations have now pushed this study well into the 
perturbative regime. Fig. [7j shows the inclusive charged hadron spectra as a function of 
centrality for Au+Au collisions at 200 GeV, extending to p T = 12 GeV. PHENIX has also 
shown 7T° spectra from both Au+Au and p+pJH Fig. || shows Raa(pt ) for these data. 
Suppression factors of 4-5 are observed in central collisions for both 7 T°S and charged 
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Fig. 8. Raa(pt) (cq. |3] or slight modification) for charged hadrons (left panel, STAR0) 
and 7r°(right panel, PHENIX0). 


hadrons, with weak, if any, dependence on p T above p T ~ 5 GeV/c. At asymptotically 
high jet Et the relative energy loss should be negligible and Raa(pt) should return to 
unity. No evidence of this limiting high energy behavior is seen. 

Fig. |9] shows a pQCD calculationil 
incorporating Cronin enhancement, nu¬ 
clear shadowing and energy loss, com¬ 
pared to the data for central collisions 
from Fig. (8j as well as lower energy 
fixed target data. The RHIC data are 
compatible with a gluon rapidity density 
dN 9 /dy ~ 800 —1200, several times that 
of cold nuclear matter. The calculation 
predicts a suppression factor independent 
of p T , an effect specific to RHIC energies 
due to the interplay between the Cronin 
effect and energy loss. In nuclear colli¬ 
sions at the SPS (fixed target), energy loss 
is masked by a much larger Cronin effect 
due to the much more steeply falling in¬ 
clusive spectrum, while at the LHC par- 
tonic energy loss overwhelms other nuclear effects. The clear separation of Cronin 



Fig. 9. pQCD calculation of hadron sup¬ 
pression, 0 compared to data from Fig. 0. 
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enhancement from energy loss at RHIC awaits the upcoming data from d+Au colli¬ 
sions. 


4.2 High pt : Elliptic Flow 


Non-central nuclear collisions generate a transversely asymmetric reaction zone. In 
high multiplicity events the orientation of the reaction plane can be determined experi¬ 
mentally through measurement of the azimuthal anisotropy of the reaction products. E! 
The occurrence of an anisotropy in momentum space (“elliptic flow”) is prima facie 
evidence of multiple interactions amongst the reaction products, which convert spatial 
into momentum-space anisotropy. Furthermore, the asymmetry is generated early in 
the collision and expansion of the system will only dilute it.il The agreement in detail 
of hydrodynamic calculations with elliptic flow measurementsil for pr < 2 GeV/c (in¬ 
cluding mass dependence, a sensitive test) argues strongly in favor of local equilibration 
early in the collision. 

Finite azimuthal correlations of high p T hadrons with the reaction plane can result 
from partonic interactions in the medium.il A priori the di-jet axis has no correlation 
with the reaction plane orientation. If the outgoing parton or its hadronic fragments 
interacts with the surrounding medium, the energy loss will depend upon the length of 
matter transversed. The consequence will be an azimuthally varying hadron suppres¬ 
sion, correlated with the reaction plane. 

The strength of this correlation is measured by the elliptic flow coefficient!! v 2 (pr) = 
(cos(2</>)), the second coefficient of the Fourier decomposition of the azimuthal distri¬ 
bution relative to the reaction plane. Fig. m left panel, shows v 2 (pr) for charged 
hadrons measured by STAR from minimum bias Au+Au collisions at 130 GeV, to¬ 
gether with a calculation patching together hydrodynamics at low p T and pQCD with 
partonic energy loss at high p T JH The hydrodynamic calculation agrees at p T <~ 2 
GeV/c. There is rough agreement with the pQCD calculation incorporating large initial 
gluon density at p T ~ 4 — 6 GeV/c. The right panel shows v 2 (pt) for the 200 GeV 
data, extending the measurement to pr ~ 12 GeV/c. The finite correlation strength per¬ 
sists for all but the most central collisions to pr ~ 10 GeV/c, well into the perturbative 
regime. Shuryak has pointed out!! that the the plateau values of v 2 for p T > 2 GeV/c in 
Fig. [TO] are in fact very large, exhausting (or even exceeding) the asymmetries expected 
from the initial geometry of the collision. 

The reaction plane for the data in Fig. [H] is calculated using hadrons with p T < 2 
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p T (GeV/c) 


Fig. 10. Elliptic flow v 2 (pt ) for charged hadrons. Left panel: minimum bias at 130 
GeV compared to pQCD calculation.ilil Right panel: centrality dependence at 200 
GeV.ll 

GeV/c. This measurement may be biased by multiparticle correlations unconnected to 
the reaction plane (“non-flow” effects), such as resonance decays, momentum conser¬ 
vation or jet fragmentation. Atp T ~ 4 GeV/c, non-flow contributes ~20% to v 2 (pr)® 
The magnitude of non-flow at higher p T remains an open problem, though preliminary 
studies indicate that the non-flow relative contribution does not increase markedly at 
higher p T . 

4.3 Highj9 T : Two-particle Correlations 

The two previous sections discussed inclusive hadron production and the correlation of 
single leading hadrons with the reaction plane determined from the bulk of the event 
(pt < 2 GeV/c). While jet fragmentation is reasonably expected to dominate hadron 
production at sufficiently high p T , it remains an open question whether the current 
measurements have achieved that limit. Full jet reconstruction in Au+Au events is not 
possible due to the complex underlying event, but intra-jet correlations amongst high 
p T hadrons are still visible and we will use these to map the transition in p T from soft to 
hard physics. The rate of back-to-back di-jets (di-hadrons) may be especially sensitive 
to partonic energy loss effects: if one escapes from the surface, its partner has enhanced 
probability to plow through the bulk matter. 
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(a) STARTS large |Ar/| 


subtracted. 


(b) PHENIxS 


Fig. 11. Azimuthal distributions of high pr hadron pairs demonstrating intra-jet corre¬ 
lations. See text for details. 


4.3.1 Near-Angle Correlations 

We first consider near-angle correlations at high p T to look for jet-like correlations. 
In the previous section, however, we showed that high pr hadrons are azimuthally 
correlated with the reaction plane (elliptic flow) with strength n 2 , so that even in the 
absence of jets, high p T hadron pairs will exhibit a correlation with strength v\. Elliptic 
flow and other non-jet contributions must be accounted for to extract quantitatively the 
correlations due to jets. 

Fig. [TT| shows the distribution of relative azimuthal angle of high p T charged 
hadrons with respect to a trigger: 

• Fig. |1 1 (a)| : the trigger is a charged hadron with p T > 4 GeV/c in the STAR 
TPC. Correlations due to non-jet effects, in particular elliptic flow, have long 
range in Ap (difference in pseudorapidity between the pair), whereas intra-jet 
correlations have short range in Arj (caveat: also true of correlations due to res¬ 
onances). The contribution of non-jet effects to the angular correlation can be 
assessed by comparing the azimuthal distributions for pairs having |A? 7 | < 0.5 
and 0.5 < | A?/| < 1.4. The figure shows the difference between azimuthal distri¬ 
butions for small and large Arp and the peak at Ao = 0 in the figure is therefore 
short range in both Ao and Arj. (The away-side peak due to di-jets is subtracted 
by construction.) Distributions for Au+Au collisions at 130 and 200 GeV are 
shown, as well as a similar analysis from UA1 for p + p at ^=540 GeV (though 
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Fig. 12. Azimuthal correlations of high p? hadron pairs in Au+Au compared to p+p 
plus elliptic flow (eq. |]).@ Left panel: peripheral collisions. Right panel: central 
collisions. 

without the large | A//| contribution subtracted). The peak has similar width in all 
cases and grows in strength with increasing y r s, as expected due to the hardening 
of the underlying jet Et spectrum with higher ^/s. 

• Fig. |1 1 (b)| : trigger is a neutral cluster in the PHENIX Electromagnetic Calorime¬ 
ter with energy > 2.5 GeV, principally photons from high p T 7r° decays. The 
azimuthal distribution from the PYTHIA event generator matches correlations 
from p+p events measured in the same detector, assumed to be due to jets, and 
PYTHIA is therefore used for comparison to Au+Au. The dashed line shows the 
expectation from elliptic flow, and the solid line shows the sum of PYTHIA and 
ellitpic flow. The near-angle correlation is clearly dominated by jet contributions. 

As an additional check, the near-angle correlation exhibits a prefered charge order¬ 
ing!^ in quantitative agreement with jet fragmentation studies at LEP.0 The quanti¬ 
tative agreement suggests that hadron prduction at px > 4 GeV/c is dominated by jet 
fragmentation in both p+p and Au+Au. 

4.3.2 Back-to-Back Hadron Correlations 

We now turn to high p T hadron pairs with large azimuthal angle difference (back-to- 
back pairs) and search for the effect of partonic energy loss on the di-jet (di-hadron) 
production ratc.ll 11 Fig. [12] shows the azimuthal angular distribution between pairs for 
the most central and peripheral Au+Au collisions at 200 GeV. The contribution of ellip¬ 
tic flow and other non-jet effects must again be accounted for. Since we are interested 
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in the correlation strength of back-to-back pairs, which have broad correlation in A rj, 
the trick used for the near-angle correlations of subtracting the large A 77 distribution 
will not work. Instead, STAR compares the full azimuthal distribution from p+p and 
Au+Au events measured in the same detector. The null hypothesis is that high p T corre¬ 
lations in Au+Au correspond to an incoherent superposition of the jet-like correlations 
measured in p+p and elliptic flow. This is expressed in the following ansatz for the 
azimuthal distribution in Au+Au: 


D 


AuAu 


= Df + B{ 1 + 2 v\ cos(2A0)). 


(4) 


Z) 2 P is the azimuthal distribution measured in p+p, v 2 is determined independently (see 
Fig. [T0|) and is taken as constant for p T > 2 GeV/c, and B is an arbitrary normalization 
to account for the combinatorial background, fit in the non-jet region A</> ~ n/2. 


< 0.75, 4<p r (trig)<6 GeV/c 
> 2.25. 4<p j1lrig)<6 GeV/c 


< 0.75. 3<p T (lrig)<4 GeV/c 
> 2.25. 3<p j(trig)<4 GeV/c 


12 



The solid lines in fig. 
show the elliptic flow compo¬ 
nent and the histograms are the 
result of eq. [4j. At all central¬ 
ities, the near-angle peak is fit 
well by eq. as is the large an¬ 
gle peak for peripheral collisions 
(but slightly less well, see be¬ 
low). Thus, back-to-back hadron 
pair production in peripheral nu¬ 
clear collisions is well described 
by an incoherent superposition 
of elliptic flow and jet-like cor¬ 
relations measured in p+p. How¬ 
ever, the large angle peak in the most central collisions is not described by eq. 0. We 
are therefore led to a striking observation: back-to-back jet production is strongly sup¬ 
pressed in the most central Au+Au collisions. 

The full centrality dependence of this effect is seen by subtracting the elliptic flow 
component and normalizing the integrated correlation peaks in Au+Au to those in p+p: 


Fig. 13. Iaa vs centrality (large N part =central 
collisions) for near-angle (upper) and back-to-back 
(lower) angular correlations.0 Left: trigger p T > 4 
GeV/c. Right: trigger^ > 3 GeV/c. 


rf(A0)[D 2 AuAu - B( 1 + 2v\ cos(2A0))] 


Iaa{ A</>i, A0 2 ) = 


(5) 


/££■ d(A4,)Df 

Fig. [13] shows I,\a as a function of event centrality for two sets of trigger and as 
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Fig. 14. Ratio of proton to pion yield vs pr for central (filled) and peripheral (open) 
collisions .0 


sociated particle thresholds. The near-angle peak strength is near or above unity for 
all centralities and both thresholds: there is no suppression in the near-angle correla¬ 
tion. In contrast, the back-to-back correlation strength in Au+Au relative to p+p de¬ 
creases smoothly from peripheral to central collisions. For the most central collisions 
the strength is consistent with zero, which would correspond to a complete suppression 
of the away-side hadron correlation rate above the threshold. 

Various nuclear effects may contribute to the features seen in fig. [13]. Initial state 
multiple scattering might generate both the increase of the near-angle and suppression 
of the back-to-back correlation strength for more central collisions. The suppression 
in the most peripheral bin for Au+Au relative to p+p for both the near-angle and back- 
to-back peaks could be due to nuclear shadowing, or possibly the interplay between 
multiple scattering and absorption in matter. These issues will be clarified by the up¬ 
coming d+Au run at RHIC. 


4.4 High p T : Baryon Enhancement 

Near-side azimuthal correlations indicate that hadron production at pt >~ 4 GeV/c 
is dominated by jet fragmentation at all centralities. However, Fig. [M] and [L5] from 
the PHENIX Collaboration showing the relative populations of mesons and baryons at 
high pt present a challenge to this picture. In fig. [14] the proton/pion ratios for periph¬ 
eral collisions are in the vicinity of 0.3-0.4, consistent with many measurements of jet 
fragmentation in elementary collisions, but the ratios for central collisions are near or 
exceed unity for p T ~ 2 — 4 GeV/c,0 very unlike conventional jet fragmentation in vac¬ 
uum. This could be due to transverse radial flow (hydrodynamic behavior), which has 
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a strong mass dependence. 1 ! extending “soft” physics into the nominally pertrubative 
regime. A second possibility is the Cronin effect: significant enhancement of baryon 
relative to meson yields for few GeV/c has been measured in fixed target p+A 
measurements! 1 ! and may be attributable to initial state multiple scattering, though this 
measurement must be repeated at RHIC energies for quantitative comparison. More 
exotic mechanisms such as baryon production via gluonic junctions^! have also been 
discussed in this context. 


PHENIX Preliminary 


The ratios in fig. [1J are measured to 
Pt ~ 4 GeV/c. It is important to estab¬ 
lish whether the anomalous baryon en¬ 
hancement for central collisions extends 
to higher pr or whether the conventional 
jet fragmentation ratio is re-established. 

As a step in that direction, fig. |T3j shows 
the ratio of 7r° to all charged hadrons for 
Pt < 9 GeV/c, for minimum bias Au+Au 
events. Even for minimum bias the ratio Rg 15 Ratio of „(> l0 total charged hadron 
is suppressed relative to the conventional yidd vs [n fro minimum bias Au+Au col . 
expectation of 0.7-0.8. This presents an dsions at 2 00 GevS 
important puzzle, but more detailed cen¬ 
trality dependence and explicit measurement of the baryon and meson fractions would 
help to ellucidate this phenomenon. 
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4.5 High Pt : Discussion 

The foregoing sections reveal the following facts about high pr hadron production in 
Au+Au collisions at RHIC energies: 

• Inclusive hadron production is suppressed in central collisions by a factor 4-5, 
with no strong p T dependence within 5 < p-p <12 GeV/c. 

• For non-central collisions, finite azimuthal correlations of leading hadrons with 
the reaction plane persist to pr ~ 10 GeV/c. 

• Near-angle two particle correlations show clear jet-like correlations for pr > 4 
GeV/c. 


19 




















• Back-to-back two particle correlations show a striking suppression in central col¬ 
lisions. 

• The baryon/meson ratio for pr < 4 GeV/c is strongly enhanced in central col¬ 
lisions, whereas it is consistent with jet fragmentation in vacuum for peripheral 
collisions. 

All but the last of these phenomena suggest a picture in which the system generated in 
nuclear collisions at RHIC is largely opaque to high energy partons and only those jets 
produced on the periphery of the reaction zone and heading outwards survive and are 
observed. This scenario naturally generates the suppression of the inclusive spectra and 
the space-mometum correlation that produces v 2 (pr), as well as the strong suppression 
of back-to-back pairs. 

Additional support for 
this picture comes from Fig. 

B which shows the yfs 
dependence of the inclusive 
yields and v 2 (p T ) between 
\/snn=130 and 200 GeV. 

The lines in the upper fig¬ 
ure show the pQCD expec¬ 
tation of the growth of in¬ 
clusive cross section with 
a/s for p+p and Au+Au 
(with and without energy 
loss). The inclusive rates 
grow by a factor 2 at p T ~ 

6 GeV/c, consistent with 
expectations from pQCD, 
whereas v 2 (pt) is indepen¬ 
dent of a/s'xn at the 5% level 
for p T > 2 GeV/c. The inclusive rates grow but the anisotropy does not change, sug¬ 
gesting that the origin of v 2 (pr) is geometric, not dynamic. 

However, this picture as presented is only qualitative, and it remains to be shown 
whether a surface emission model can simultaneously describe the inclusive suppres¬ 
sion, v 2 (pt), and the suppression of back-to-back pairs. The anomalous baryon/meson 
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ratio presents a puzzle and does not seem to fit into this picture. That effect may be 
due to conventional mechanisms such as transverse radial flow or initial state multiple 
scattering, or more exotic physics such as baryon production via gluonic junctions.El 
Further analysis of existing data and measurements of d+Au collisons at RHIC are 
needed to clarify these issues. 

Finally, it is necessary to establish what, precisely, is interacting with the medium 
and being absorbed: partons or hadrons? Simple estimates of hadron formation time 
give tf ~ ( E/m)r for a hadron of radius r,@ so that a pion with p-r ~ 5 GeV/c is 
formed well outside of the hot and dense medium. Alternatively, string interactions in 
matter prior to fragmentation lead to suppression the leading particle yields and a much 
shorter apparent formation time for high pt hadrons.il 

5 Summary and Outlook 

RHIC has had two successful data-taking periods. It has met its design luminosity goal 
for Au+Au, and has operated as the world’s first polarized proton collider. Upcoming 
runs in 2003 will study asymmetric (d+Au) collisions for the first time in a collider, 
together with an extended polarized proton run. A long Au+Au run at top energy will 
occur in 2004, possibly with an energy scan at lower integrated luminosity. 

An extensive array of measurements in the soft physics sector indicate the pro¬ 
duction of an equilibrated system at high temperature and pressure which blows apart 
rapidly. In the high p T sector, striking signals have been observed that suggest strong 
partonic energy loss resulting from a system at high energy density at early time, though 
there remain important puzzles and open questions. 

I have not discussed the study of heavy quark production, which is in its infancy at 
RHIC. Charmonium suppression, a promising signature of deconfinement,i has been 
observed in fixed-target nuclear collisons at the SPS.0 PHENIX has measured the 
inclusive electron spectrum, due primarily to charm,0 and recently reported the first 
measurement of J/-0 at RHIC.E3 This physics will be a main focus of future runs. 
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